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ABSTRACT 

This paper presents approximate analytic expressions for the emissivity 
and absorption coefficient of synchrotron radiation of mildly relativistic 
particles with an arbitrary energy spectrum and pitch angle distribution. From 
these, an expression for the degree of polarization is derived. The analytic 
results are compared with numerical results for both thermal and non-thermal 
(power law) distributions of particles. 
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I. INTRODUCTION 


The formulas for evaluation of emissivity and absorption coefficient of 
synchrotron radiation in the ultra-relativistic (synchrotron) and non- 
relativistic (cyclotron) limits have been known for decades.^ In the 
intermediate energy range, however, no simple formula exists for an arbitrary 
distribution of particles. Earlier Trubnikov derived formulas for emissivity 
and absorption limited to a thermal gas and convenient only for propagation 
perpendicular to the magnetic field. As a result, the usual practice has been 
to use lengthy numerical calculations.'^* . 

In a recent paper^ (referred to as Paper I, hereafter) we presented 
simple approximate methods for evaluation of the frequency spectrum and 
angular variation of the synchrotron radiation at high harmonics from an 
(essentially) arbitrary distribution of particles in a given magnetic field, t 
this paper, we shall use the same methods to calculate the emissivity and 
absorption coefficient of the extraordinary and ordinary modes of synchrotron 
radiation separately. Also, we will derive the degree of circular polarization 

• i 

from these expressions. The general expressions presented in Sec. II were 
derived for high harmonics of the cyclotron frequency which are the most 
useful since the low harmonics are usually self-absorbed, absorbed by the 
surrounding plasma, or suppressed by the Razin-Tsytovich effect. In Sec. Ill 
we use the above results to find emissivity, absorption and polarization of 
the radiation from particles with Maxwellian and power law energy 
distributions. Here, through comparison with numerical results, we also show 
that our formulas provide a good approximation even at low harmonics. In Sec. 
IV we present a final summary. 
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II. GENERAL RESULTS 
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Consider particles with charge e, mass and distribution £(U»Y) where 
fdpdy is the nuabev density of particles in the energy interval (in units of 
m^c^} from y to y -t- dy and with pitch angle cosine between p and y + dp. 
The emissivity and absorption coefficients for the ordinary {*) and 
extraordinary (-*) modes at frequency V and at angle 6 with respect to the 
magnetic field B are^*^. 



where 

1 *“ ( X J^(x .) 

o — _ — y ) T (cos0-3p)j (x ) - (1 - PpcosO) 

m=l ( “ 

3cos0-y 

6) *• “ 31ntf(p,y)/PY*3/9Y + T7:: ainf(p,Y)/3p * (3) 

&Y*(1“P*)^'^ 

y ® - (1 - ppcos0), X. “ (vyv^)PsinO (1 Vj, ® eB/2m^c . (A) 

The quantity T represents the ratio o£ the major to the .minor axis of 
the polarization ellipse and, in general, is a complicated function of angle 
0, V, and the plasma frequency Vp. The ordinary (+) and extraordinary 

(-) modes are distinguished by their respective values of T, which in vacuum 
or for Vp/v « 1 are 

T^ - ~T_^^ “ ^ K " V^sin^e/CZucose) . (5) 
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Note that for >> 1, the quantity ^ << 1 and T+ ±1 except for a 

small range of angles ( it/ 2 - 9) < Vj^/2 where C > 1. And in the limit 0 **• 
7 t/ 2, C “ ®o that T^*>- » and T_-^0. 

Since we are only dealing with the higher harmonics, m is large and ve 
can replace the sum in equation (2) with an integral. This integral can then 
be done using the 6- function. Also, since m is large, ve can use the 
Carlini approximation for the Bessel's function as^. 


J (mz) 
m 






^2a-z2)'* 


( 6 ) 


which is valid for « 1. 

Ve evaluate the integrals to first order in this quantity. 

Integration over Pitch Angle 

As can be seen from £q. (6), the Bessel function and its derivative hav^i' 
sharp maxima at y ° 3cos9 and tend to zero rapidly at y^^^l. The rest of 
the integrand is positive and var. ''s slowly with y if the pitch angle 
distribution of the particles is not extremely anisotropic. We use the method 
of steepest descent to evaluate the integrals over y. The only exception to 
this is the coefficient of iii Eq. (3) which becomes zero at y 

“ COS0/3. For a general cos0 this will give rise to a secondary maximum whose 
contribution is of second order in our expansion parameter. In particular, 
for perpendicular propagation (cos0* 0) this term becomes zero at the 
maximum of the function Z*“. This gives rise to two maxima, each with a second 
order contribution to the integral. 

The use of the method of the steepest descent amounts to setting 
y = 3cos9 in the integrand and multiplying it by - 3^cos*0 )3/^. 
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This means that most of the contribution to the radiation at an angle 6 

comes from particles with pitch angle cosines in a narrow range of the 

order of (frv centered at y " 3cos9. 

b 

This approximation is valid only for pitch angle distributions which are 
not extremp\y anisotropic, i.e., distributions with 2lnf/3y << vA>jj Icf. 
Paper 1, Eqs. (7) and (19)1. Note that this enables us to drop the second 
term in the expression for (o in Eq, (3) so that the pitch angle dependence 
becomes the same for j and ic. 

Integration over Energy 

The resultant integrand for integration over energy also has a sharp 

maximum for particle distributions which fall rapidly with increasing energy 

"6 

[e.g. , a power law f « (y - 1) or a thermal distribution f « exp(-y/kT)]. 
Consequently, we can use the method of steepest descent to carry out this 
integration also» This gives 


u 


2 


r 

ite / 



c [y 

'J 1 


2vu /v. Y 
' o b'o 






X (y )f (3 cos0,Y ) 
K *0 K o 'o 


(7) 


where 




Y(e,y) 


Vd^lnf^/dy^ - 



WYdlnfj../dy 

“ 1 + P^Y^sin^O 

y (1 + T^) . 


( 8 ) 

(9) 

( 10 ) 
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Mote that all these expressions are evaluated at critical energy (and 
the corresponding and “ 1 + sin^6), where most of the 

contribution to the integral comes from. There are two such critical 
energies; one for emission and one for absorption. These are obtained from 
the transcendental equation; 


original page fS 

OF POOR QUALnV 


u"^(u^ - 1)“* + (1 - cot^e/yOlnZ^ - -0. /2vsin’0)41n£. /dv 

D J »K 

= e/sin* 6 . 


( 11 ) 


0 

Here and in equation (7) the function f • and f are related to the 

J ^ 

particle distribution f as 


f. “ f(0cos6,Y)/Y . f - -Mlf(6cos0,Y)/3Y^l/<iY • (12) 

J 


In Eq. (8) the function H is a complicated function of and 0. 

Hovever> as shown in paper 1, W can be approximated by the following simple 
expression, 

W - 3/2 + 1/(y^ - 1), (13) 

0 

which is an excellent approximation in the entire region where equations (6) 
to (11) are valid. 

Ordinary and Extraordinary Modes 

The emission and absorption coefficient of each mode is obtained by the 
substitution of T^ in Eq. (7). Thus, and become proportional to 
Y+ evaluated at the respective values of the critical energy y^. Since T^T_ « 
- 1 , it can be shown that 

y+ " (u^-2T^ucosO + T^^cos^6)/lY^sin* 0 (1 + T^^)J, 

Y- “ 2T^ucos 0 + cos*0)/lY*sin*0(l + T^^)J. 

6 


( 14 ) 
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The total emissivity then becomes 


itot “ i+ ♦ i- ® cos*0)/Y“slne - 1+2 cot*e/Y* (15) 


in agreement with Paper I. Similarly, with the help of Eq. (5) the difference 
of emisaivities is given by 

- (2ucoae/Y*sin 0 + 0/(1 + (16) 

For V /v. » 1 and away from the direction perpendicular to the magnetir. 
b 

field, L 1 ^ ^ that the ratio of j and K of the ordinary 

to that of the extraordinary mode becomes 


K+/»c„ - j+/j_ 



/ + ^os2q 

\ - cos^e / 


1 - 2C 


(17) 


This approximation breaks down for (0 - tt/2) <Yq/2v. As 0 -*■ ir/2, u Y » 

Tj, = 2.C -► « and 


(u - Vjj/v)^ V*COS*0/v^Y* 
Y_ 1 + 2vcot®0/^)^Y 


(18) 


In the limit 0 » ir/2, « 0 and Y_ “ 1. This, of course, means that for 

evaluating the and at 0 " tt/2 one must consider higher order terms 

which we have neglected. At 0 * ir/2 the ratio Y^/Y_ is of higher order of 
our expansion parameter n*'*(l“Z^)"” ' “ ^oT'o^b^'^' 
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1 . Polaritation 

The polarization of the radiation can be obtained from the Stokes 
parameters which are (cf. references 1 and 4) 

Jq - + «cos6T^(jJ_)l/2l(l + tJ) . 

iu " 2sin«(j^j_)^^^ , (19) 

With the help of Eq. (5) and using Eqs. (15) and (16) and if we define 

Y^-Y„ (2uco8O/Y*sic.*0+5) 

p m — — w — — — — ( 20 ) 

^ Y^+Y_ (l+U^^^(l+2cot^e/Y*) 

these become 

V-jtot “ 1^1*0 ^ . 

iu^itot " ■* 

jy/itot - 1^0 - • 

Here £ is the phase difference between the ordinary and the 
extraordinary modes* For a definite value of 5 the degree of Polarization P 
» (jq^ + ^ 0 ''®ver, for radiation from many particles 

in a large source with large Faiaday rotation* the phase relations are 
randomized so that the average values <sin6> " <cos5> “ 0. In this case jy ■* 
0 and the degrees of linear* circular and total (elliptical) polarization are 
simply 


Qualhy 


•■u* - 'circ ■ V“ * '■tot “ ■’o- 


( 22 ) 


8 



When v/Vjj » 1 and ^ « 1, « ^circ radiation is circularly 

polarized except near 6 " ir/2 where and radiation becomes linearly 

polarized. As evident from Eqs. (20) and (22) the degree of polarization 
decreases with increasing frequency (or the energy of the particles). 

2. Sumittarv 

Eqs. (7) to (13) along with (21) and (22) give our results in their 
general forms and are valid for all particle distributions which are not 
extremely anisotropic. The results for the extremely anisotropic situation are 
more complicated and were described in Paper 1. The modification for 
obtaining the emission and absorption coefficients of the two modes separately 
and the polarization is similar to the modification described above. We will 
not present these results here because of their complexity and their limited 
usefulness. 

Given a distribution function subject to this limitation, the first step 
is evaluaticn of the critical energies from Eq. (11). Then Eqs. (6) to (10) 

and (12) and (13) evaluated at the appropriate y^s give the desired 
results. The most complicated part of this procedure is the solution of Eq. 

(11) for Yq • turns out that for most practical cases it is not necessary 
to solve this equation. 

In the next section we shall show how this step of the calculation is 
simplified considerably for the two most commonly used particle distributions. 
Before doing so we consider the asymptotic limits of these equations. 

3. Asymptotic limits. Let us consider first the case when angle 6 is not 
too small (i.e., radiation away from the direction of the field). Then in the 

two extreme limiting cases, Eq. (11) simplifies t '- * ’ * 

OF POOK 

i) e « 1, Y 3*Y* “ 2/3esinG , W " 3/2; 

' o o o 

(23) 

ii) E » 1, P^«l» 1, PqYo '•A/e, w« e/4. 
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The first case is realized at high frequencies and particle distributions 
which are not extremely non-relativistic which is the case of interest here. 
The second case is valid for non-relativistic particles and at low frequencies 
and has limited usefulness except for a low temperature thermal gas. 

III. EMISSIVITY AND ABSORPTION COEFFICIENT OF TWO COMMONLY USED 

PARTICLE DISTRIBUTIONS 

The two particle distributions we use as examples are i) the distribution 
from a thermal gas, i.e. a Maxwellian distribution in energy and isotropic 
pitch angle distribution; and ii) the distribution with a power law spectrum 
at high energies and with a slowly varying pitch angle distribution. 

A. Thermal Spectrum 

In this case, the distribution f of particles at temperature kT (in 
units of BgC^) is 

£((.,Y) - Ce-(T-l)/kIY ( Y* . 1 )i/2 . (24) 

where for kT < 1 

C (n/2) l 2 Tr(kT) 3 j-i /2 q _ 15 + ...). (25) 

and n is the number density of particles. From these and Eq. (12), it is 

clear that kT^ “ fj ■ f/y so that [as is evident from Eq. (11)1 the 

critical y is the same for both j and k and, in fact, according to Eq. 

0 

(7), “ kT<^. We also find that 

d1n(f/y) -1 / kTy \ d^ln(f/y) -(y* + 1) 

_ [1 1 ^ , (26) 

dy kT \ yM/ dy* (y^ - 1)^ 
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Using these equations we can calculate Yq X from Eqs. (11) and (8) 

respective!). As shown in Paper these expressions can be considerably 

simplified. We find that the following expressions . ta 
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. j (2vkT/\^)Cl + 4.5 vkTsin%/Vjj)“^^^ kT < 1 



(4VkT/3Vj,sinO)^/3 


kT 1 


X 2 - (2kT/Y^)( Y* - - 1), 


kT <: 1 

'U 


(27) 


I 


have the correct asymptotic limits in agreement with Eq. (23) and agree with 
the exact results from Eqs. (8) and (11) to within 30Z for most relevant 
ranges of anglest frequencies and temperatures and better than lOZ in the 
majority of the interesting cases. 

The above equations and Eqs. (7)» (8), and (10) give a complete 
description of the emissivity and absorption coefficient from a Maxwellian gas 
at all temperatures and frequencies. They are valid for kT .v 1 because at 
temperatures kT > 1 the use of the method of steepest descent for integration 
over the energy becomes less accurate. However, the existence of the 

O 

extremely relativistic thermal gas is in doubt . On Fig. 1 we compare the 
total absorption coefficient k “ <+ + k _ obtained from these relationships 
with numerical results from Lamb and Masters'^. As evident, our analytic 
results give excellent agreement t> the detailed numerical results even at low 
harmonics. A more detailed comoarison with similar results was presented by 
Marsh and Dulk^. 

In the two limiting cases described in the previous section, these 
equations are considerably simpler. The interesting case, e << 1 corresponds 
to V kT/Vjj » 1, Yq “ 4vkT/3VjjSinO, so that 


11 


( 28 ) 


jj. Atk'j- .xpLl ---- 


Y . 1 ♦ 2(3V .in‘e/4vkT)l/3 OniGfNAL PAGE [9 
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Similarly^ the degree of polarization becomes P - (48Vj,ain'0/vkT)^^^ 


B. Power Law Energy Spectrum 

Power lav spectra are commonly used spectra in astrophysical problems and 
in other problems when the tail of the Maxwellian distribution begins to 
deviate from the exponential form. Usually power law spectra are defined with 
a low eutergy cut-off. To avoid such discontinuities and the divergence of the 
number of particles, we assume a spectrum of the form 

n(6 - i) * 

£(u y5 “ — t:: W “ i)/ c )f g(u)» 

Here plays the role of the low energy cutoff (in units of mgC ). For 

energies much greater than the spectrum is a power law with index ~c 

but it tends to a constant value at lower energies. The particles can be 
classified ' as ultra-relativistic or non-relativistic if e » 1 or c « 1. 
We are interested primarily in cases with c^=l. 

For distributions which are not highly anisotropic (i.e. dlng(p)/dp << 
v/V|j}, we can carry out a calculation similar to that for a thermal gas. From 
Eqs. (12) and (29) we find that 


/*.« 


g(y)dp - 1 


6 2y*- 1 

Y- 1 Y(Y^- 

so that the quantities and X ^ in Eqs. (8) and (11) can all be calculated 
for exact evaluation of the emissivity and the absorption coefficient in Eq. 
(7). 

As shown in Paper I for semi-relativistic particle energies, 1, and 

for vA> >>1 all these complicated expressions can be simplified 
b 
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cc ixsiderably. This is because it turns out that the high energy asymptotic 
limit of Eq. (23) provides a good approximation throughout most of the 
relevant ranges of frequency and angle 6. In this licit Eqs. t29) and (30) 
give 


PV = X^4v/(3v^sin9), » 

0 0 D 


1+6 for j 
2+6 for K 


(31) 


Values of obtained from this simple expression agree to within. 30% of 

the exact values derived from Eq, (11) for 1 ^ ^ 0.5 and sinO > , Note 

also that in this limit the exact form of the energy distribution is not 

important as long as it tends to a power law f « y high y. In the 

extreme-relativistic limit v >> v, and Y >> 1> substitution of Eq. 

bo ' 

(31) into Eq, (7) gives 


j±(v.O) 

(v,6) 




•s/Sire V, sinO 
b _ _ 

Ac 


C D 


Ce g(3 cosO)Y (0) 
CO - 


(32) 


(3£*v (6+2)sin0/4v) e 
c b 


6/2 -( 6 + 2)/2 


Pf)- 


Note that in this limit Y 4 . “ 1 + 2 cot9/y so that the ratios 

Z. o 


of 


«L/k:^ = “ O + ^ cot0/y^) 


(33) 


tend to unity with increasing frequency. This, and the dependence on V/Vj^, 
:i»id 6 of Eq. (32) is identical to the results for emissivity obtained from 
ultra-relativistic expressions. 

It is interesting that with a quite different method and approximation 

we have obtained the same expression. The reason for this agreement can be 

I “3/2 

seen by examination of the expansion parameter (1-z) /m= uyv /v. 

00 b 

which in this limit is equ-.l to 4/13(1 +6)J or 4/[3(2 +6)J. For 6 > 3 the 
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expansion parameter is less than 0.3 which, although not extremely small, 
nevertheless is less than unity. Thus, for an ultra-relativistic expression 
more accurate than 30Z to 502, one must include higher order terms in our 
analysis. 

For extremely non-relativistic particles, that is, for « 1, the above 

expressions are valid as long as \f^6/ V « 1. For the unlikely case of Vj^6/V 

>> 1, we find tnat ~ ^ ~ Substitution of this in (7) 

gives the emissivity identical to that expected from a thermal gas if one 

identifies t^/6 with the temperature kT (cf. Paper 1). 

In Figs. 2, 3 and 4 we compare results obtained from substitution of Eqs. 

(30) and (32) in Eqs. (7) and (20) for j^, , and P, at 0 60° and 

5 ° 4 with the results of numerical integration kindly provided by Marsh and 

8 

Dulk , who used a power law spectrum with a sharp cutoff at 1.02. 

The values for and (c^ are within 60% at high frequencies and the 

approximation is better at higher v/v^^. The polarization looks better at 
lower frequencies (although it still is better than 50% even when P is small 
and errors can be magnified). Also our results are (nearly) systematically 
higher than the numerical ^su ts at high harmonics, which Could be the effect 
of higher order terms. However, the percentage error never exceeds 30% for 

v/Vjj ^ 6 . 

IV. SUMMARY 

Using a simple method of integration developed previously, we have 
derived expressions for the emissivity, absorption coefficient and 
polarization of synchrotron radiation for an arbitrary distribution of 
particles. 

Equations (7) to (13) and (21) give our results in their most general 
form. And we find that Eqs. (8) and (11) can be simplified considerably, as 
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in Eqs. (27) for a Maxwellian distribution and in (31) for a power law 
distribution. 

Our results do agree with previous analytic results, and they give good 
approximations to detailed numerical results. Although our results were 
derived for high harmonics, they give good agreement down to lower harmonics: 
to V - 6Vjj for j^, to V = l(h>b for <+» sod to even lower harmonics , V 
“ 2Vjj, for the total j and K . These results are limited to pitch angle 
distributions which are not extremely anisotropic and energy spectra that 
decrease rapidly with increasing energy. They also are only applicable for 
emissivities and absorption coefficients away from the direction of the 
magnetic field lines. 

Our equations are intended for semi-relativistic particles, but they also 
give excellent approximations for the extreme non-relat ivistic and ultra- 
relativistic particle distributions. 
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FIGURE CAPTIONS 


Figure 1. The total synchrotron absorption for a thermal source at 0 “ tt/ 2 
for kT ■■ 0.04 (20 keV electrons). Points are from analytic expressions; the 
solid lines are numerical results of Lamb and Masters3 

• 

Figure 2. Synchrotron emissivity of each mode divided by magnetic field B and 
total particle number N. Log (j^./BN) vs Log (v/Vj,), at 6“ 60°, 6 » 4. 

( in units erg(cm^ sec sterad H The ordinary mode has been 

shifted down by a factor of 10 for clarity. The solid lines are numerical 
results of Bulk and Marsh (private communication). The O's are our analytic 
results. 

Figure 3. Same as Fig. 2 except for the absorption coefficient v^ic^ in units 
of erg sec ■ (cm** sterad) 

Figure 4. Degree of Circular Polarization vs Log (v/Vj^) in the limit 
C 0. Circles from Eq. (22), Solid line from numerical results. 


17 


OF POOR QVJALrr^. 



Figure 1 
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